Dynamic time-step estimates for two-dimensional transient field problems using square elements

1998 ◽  
Vol 42 (1) ◽  
pp. 1-14 ◽  
Author(s):  
Rabi H. Mohtar ◽  
Larry J. Segerlind
Keyword(s):  
Two Dimensional ◽  
Transient Field ◽  
Time Step ◽  
Dynamic Time

DYNAMIC TIME STEP ESTIMATES FOR ONE-DIMENSIONAL LINEAR TRANSIENT FIELD PROBLEMS

1999 ◽  
Vol 42 (5) ◽  
pp. 1477-1484 ◽  
Author(s):  
R. H. Mohtar ◽  
L. J. Segerlind
Keyword(s):  
Transient Field ◽  
Time Step ◽  
One Dimensional ◽  
Dynamic Time

scholarly journals A two-dimensional glacier–fjord coupled model applied to estimate submarine melt rates and front position changes of Hansbreen, Svalbard

2018 ◽  
Vol 64 (247) ◽  
pp. 745-758 ◽  
Author(s):  
E. DE ANDRÉS ◽  
J. OTERO ◽  
F. NAVARRO ◽  
A. PROMIŃSKA ◽  
J. LAPAZARAN ◽  
...  
Keyword(s):  
Coupled Model ◽  
Small Scale ◽  
Two Dimensional ◽  
Time Step ◽  
Software Packages ◽  
Front Position ◽  
Circulation Patterns ◽  
Glacier Front ◽  
Order Of Magnitude ◽  
Frontal Ablation

ABSTRACTWe have developed a two-dimensional coupled glacier–fjord model, which runs automatically using Elmer/Ice and MITgcm software packages, to investigate the magnitude of submarine melting along a vertical glacier front and its potential influence on glacier calving and front position changes. We apply this model to simulate the Hansbreen glacier–Hansbukta proglacial–fjord system, Southwestern Svalbard, during the summer of 2010. The limited size of this system allows us to resolve some of the small-scale processes occurring at the ice–ocean interface in the fjord model, using a 0.5 s time step and a 1 m grid resolution near the glacier front. We use a rich set of field data spanning the period April–August 2010 to constrain, calibrate and validate the model. We adjust circulation patterns in the fjord by tuning subglacial discharge inputs that best match observed temperature while maintaining a compromise with observed salinity, suggesting a convectively driven circulation in Hansbukta. The results of our model simulations suggest that both submarine melting and crevasse hydrofracturing exert important controls on seasonal frontal ablation, with submarine melting alone not being sufficient for reproducing the observed patterns of seasonal retreat. Both submarine melt and calving rates accumulated along the entire simulation period are of the same order of magnitude, ~100 m. The model results also indicate that changes in submarine melting lag meltwater production by 4–5 weeks, which suggests that it may take up to a month for meltwater to traverse the englacial and subglacial drainage network.

Enhanced Explicit Scheme to Solve Transient Heat Conduction Problem

2007 ◽  
Vol 2 (1) ◽  
Author(s):  
Ganesh Hegde ◽  
Madhu Gattumane
Keyword(s):  
Heat Conduction ◽  
Correction Factor ◽  
Truncation Error ◽  
Transient Heat Conduction ◽  
Explicit Scheme ◽  
Stability And Convergence ◽  
Two Dimensional ◽  
Time Step ◽  
One Dimensional ◽  
Partial Derivatives

Improvement in accuracy without sacrificing stability and convergence of the solution to unsteady diffusion heat transfer problems by computational method of enhanced explicit scheme (EES), has been achieved and demonstrated, through transient one dimensional and two dimensional heat conduction. The truncation error induced in the explicit scheme using finite difference technique is eliminated by optimization of partial derivatives in the Taylor series expansion, by application of interface theory developed by the authors. This theory, in its simple terms gives the optimum values to the decision vectors in a redundant linear equation. The time derivatives and the spatial partial derivatives in the transient heat conduction, take the values depending on the time step chosen and grid size assumed. The time correction factor and the space correction factor defined by step sizes govern the accuracy, stability and convergence of EES. The comparison of the results of EES with analytical results, show decreased error as compared to the result of explicit scheme. The paper has an objective of reducing error in the explicit scheme by elimination of truncation error introduced by neglecting the higher order terms in the expansion of the governing function. As the pilot examples of the exercise, the implementation is aimed at solving one-dimensional and two-dimensional problems of transient heat conduction and compared with the results cited in the referred literature.

scholarly journals Numerical Analysis of Unsteady Laminar Flow past a Pentagonal Obstacle

2021 ◽  
Vol 10 (2) ◽  
pp. 968-974
Keyword(s):  
Laminar Flow ◽  
Flow Direction ◽  
Rectangular Channel ◽  
Time Frame ◽  
Two Dimensional ◽  
Fluid Inertia ◽  
Time Step ◽  
The Face ◽  
Computational Fluid Dynamics Cfd ◽  
Current Article

The current article dispenses the numerical investigation of a two dimensional unsteady laminar flow of incompressible fluid passing a regular pentagonal obstacle in an open rectangular channel. The centre of attention of this work is the comparison of drag coefficients estimated for two distinct cases based on the orientation of face and corner of an obstacle against the flow direction. The numerical results shows that the corner – oriented obstacle bring about 42% larger value of drag coefficient at Re = 500 than face – oriented obstacle. The substantial growth in the expanse of vortex behind obstacle (presented as a function of fluid inertia 25 < Re < 500) is analyzed through the contours and streamline patterns of velocity field. The two eddies in the downstream become entirely unsymmetrical at Re = 500 for both the cases, whereas; the flow separation phenomena occurs a bit earlier in the face – oriented case at Re = 250. Two dimensional Pressure – Based – Segregated solver is employed to model the governing equations written in velocity and pressure fields. The numerical simulations of unsteady flow are presented for 50 seconds time frame with time step 0.01 by using one of the best available commercial based Computational Fluid Dynamics (CFD) software, ANSYS 15.0.

Fully discrete numerical schemes of a data assimilation algorithm: uniform-in-time error estimates

2019 ◽  
Vol 40 (4) ◽  
pp. 2584-2625 ◽  
Author(s):  
Hussain A Ibdah ◽  
Cecilia F Mondaini ◽  
Edriss S Titi
Keyword(s):  
Data Assimilation ◽  
Error Estimates ◽  
Stokes Equations ◽  
Time Discretization ◽  
Dissipative Systems ◽  
Complete Analysis ◽  
Two Dimensional ◽  
Time Error ◽  
Time Step ◽  
Fully Discrete

Abstract Our aim is to approximate a reference velocity field solving the two-dimensional Navier–Stokes equations (NSE) in the absence of its initial condition by utilizing spatially discrete measurements of that field, available at a coarse scale, and continuous in time. The approximation is obtained via numerically discretizing a downscaling data assimilation algorithm. Time discretization is based on semiimplicit and fully implicit Euler schemes, while spatial discretization (which can be done at an arbitrary scale regardless of the spatial resolution of the measurements) is based on a spectral Galerkin method. The two fully discrete algorithms are shown to be unconditionally stable, with respect to the size of the time step, the number of time steps and the number of Galerkin modes. Moreover, explicit, uniform-in-time error estimates between the approximation and the reference solution are obtained, in both the $L^2$ and $H^1$ norms. Notably, the two-dimensional NSE, subject to the no-slip Dirichlet or periodic boundary conditions, are used in this work as a paradigm. The complete analysis that is presented here can be extended to other two- and three-dimensional dissipative systems under the assumption of global existence and uniqueness.

A finite‐element solution for the transient electromagnetic response of an arbitrary two‐dimensional resistivity distribution

Geophysics ◽  
1986 ◽  
Vol 51 (7) ◽  
pp. 1450-1461 ◽  
Author(s):  
Y. Goldman ◽  
C. Hubans ◽  
S. Nicoletis ◽  
S. Spitz
Keyword(s):  
Finite Element ◽  
Sparse Matrix ◽  
Equation System ◽  
Electromagnetic Response ◽  
Implicit Time ◽  
Two Dimensional ◽  
Time Step ◽  
Transient Electromagnetic ◽  
Exact Boundary ◽  
The Time Domain

We present a numerical method for solving Maxwell’s equations in the case of an arbitrary two‐dimensional resistivity distribution excited by an infinite current line. The electric field is computed directly in the time domain. The computations are carried out in the lower half‐space only because exact boundary conditions are used on the free surface. The algorithm follows the finite‐element approach, which leads (after space discretization) to an equation system with a sparse matrix. Time stepping is done with an implicit time scheme. At each time step, the solution of the equation system is provided by the fast system ICCG(0). The resulting algorithm produces good results even when large resistivity contrasts are involved. We present a test of the algorithm’s performance in the case of a homogeneous earth. With a reasonable grid, the relative error with respect to the analytical solution does not exceed 1 percent, even 2 s after the source is turned off.

Transient Two-Dimensional Problems Utilizing The Least Squares Algorithm

1975 ◽  
Vol 97 (3) ◽  
pp. 467-469 ◽  
Author(s):  
J. C. Bruch ◽  
R. W. Lewis
Keyword(s):  
Least Squares ◽  
Finite Difference Methods ◽  
Computer Time ◽  
Two Dimensional ◽  
Time Step ◽  
Time Stepping ◽  
Large Time Step ◽  
Least Squares Algorithm ◽  
Methods Analytical ◽  
And Storage

The least squares time-stepping algorithm, which has previously been shown by the authors to be competitive for one-dimensional problems, is applied to the solution of several two-dimensional examples having constant material properties. The results are compared against answers obtained using recurrence relationships based on the finite element and finite difference methods. Analytical results for one of the examples are also used for comparison. The least squares algorithm proved to be more accurate for equal values of time step especially in the large time step cases. It, however, requires more computer time and storage than the other methods used. Several other limitations of the scheme are also presented.

Sign in / Sign up

Export Citation Format

Share Document